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I. ABSTRACT 

Sensors to remotely measure the lin- 
ear polarization of ground scenes have 
been proposed for the Multispectral Re- 
source Sampler (MRS) , a satellite sensor 
system proposed to compliment the Thematic 
Mapper. At present justification for a 
sensor on MRS to measure scene polariza- 
tion is limited. This paper discusses a 
model for the amount of linearly polarized 
light reflected by the shiny leaves of 
such crops as wheat, corn, and sorghum. 

The theory demonstrates that, potentially, 
measurements of the linearly polarized 
light from a crop canopy may be used as an 
additional feature to discriminate between 
crops. Examination of the model suggests 
that, potentially, satellite polarization 
measurements may be used to monitor crop 
development stage, leaf water content, 
leaf area index, hail damage, and certain 
plant diseases. The model adds to our 
understanding of the potential information 
content of scene polarization measurements 
acquired by future satellite sensor sys- 
tems such as MRS. 


II. INTRODUCTION 

A series of investigations has shown 
that remotely sensed Landsat satellite 
multispectral scanner (MSS) data can be 
used to accurately identify and measure 
hectarages of crops over large areas. Ex- 
perimental global wheat production fore- 
casts have been made by melding crop area 
estimates derived from Landsat data and 
estimates of crop yield made by regression 
models based on historical weather and 
yield data 1 ' 2 . Despite these successes, 
there are limits to the present technology: 
for example, during the Large Area Crop 
Inventory Experiment (LACIE) , there was a 
tendency for spring wheat to be confused 
with other small grains such as barley. 

And potentially in the multitemporal spec- 
tral responses of various crops there is 
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information related to crop vigor and 
growth stage, information needed for var- 
ious models to predict crop grain yield 
There is a continuing need to better dis- 
criminate between crops and a need to re- 
search and develop remote sensing data 
analysis techniques to predict crop grain 
yield 6 . 

Sensors to remotely measure the linear 
polarization of ground scenes have been 
proposed for the Multispectral Resource 
Sampler (MRS) , a satellite sensor system 
proposed to compliment the Thematic Map- 
per". If polarization sensors are to be 
included in the MRS, it should be estab- 
lished that polarization measurements of a 
scene provide nonredundant information ex- 
ceeding that already found in the scene 
spectral response, now routinely measured 
by Landsat MSS with four wavelength bands 
and soon to be measured by the Thematic 
Mapper with seven wavelength bands. 

At present justification for a sensor 
on MRS to measure scene polarization is 
limited. Using data obtained in the labo- 
ratory and with an aircraft, Egan, Egan 
and Hallock, and Egan, et al. found evi- 
dence that the degree of linear polarization 


Figure 1. Specular Reflection. The 
camera received specularly reflected sun- 
light from the bright areas tf these 
wheat flag leaves. 
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(a) (b) 

Figure 2. Polarized Light from Can- 
opy. These photos, taken with a polarizer 
oriented for transmission of maximum spec- 
ularly reflected light (a) and minimum 
specularly reflected light (b) , demon- 
strate that the specularly reflected light 
is polarized. 

of the response of a scene provides addi- 
tional discriminatory information with 
which to classify the sceneW 0 ' 11 . Egan 
reached a potentially important conclusion 
that drying of leaves generally increases 
their depolarizing properties 1 . Curran 
used a photographic measurement technique 
to relate soil surface moisture to the 
proportion of polarized light in the scene 
response 12 ' 13 . In an appendix Curran pre- 
sented data showing a possible link between 
the percent linear polarization of a can- 
opy an^ its roughness 12 . 

This paper discusses a model for the 
amount of linearly polarized light re- 
flected by a plant canopy. The mod' 1 is 
based upon the morphological and philo- 
logical characteristics of the canopy and 
upon the Fresnel equations which describe 
the light reflection process at the smooth 
boundary separating two dielectrics 1 ‘*. The 
theory demonstrates that the linear polar- 
ization response of the plant canopy po- 
tentially contains information to help 
discriminate between crops. The theory 
relates the response to the agronomic con- 
dition of the crop--to factors such as 
growth stage, plant vigor, and leaf area 
index (LAI ) . 


III. POLARIZATION MODEL 

As may be verified with a polarizing 
filter and camera, the shiny leaves on a 
plant (Fig. 1) provide the basis for un- 
derstanding how light is specularly re- 
flected and polarized (Fig. 2) by a 
healthy, vigorous plant stand. The shiny 
leaves of many plants, including wheat, 
corn, and sorghum varieties, have a leaf 
skin or cuticle covered by a wax layer 
which specularly reflects light in accor- 
dance with the Fresnel equations. Because 


the light is specularly reflected, the 
equations show that it is polarized for 
all but two angles of incidence (0 and 90 
degrees) . Breece and Holmes found that 
the bidirectional light scattering charac- 
teristics of wheat, corn, and soybean 
leaves are intermediate to the scattering 
characteristics of diffuse and specular 
reflectors, indicating that the specular 
portion is an important part of the total 
leaf response 15 . 

From the laws of optics a leaf seg- 
ment with a tendency toward specularity 
portends that a significant portion of the 
surface area is flat and similarly orient- 
ed. Indeed, there is electron micrograph 
evidence that the wax deposits on glossy 
leaves may form smooth films on the cuti- 
cle or platelets which lie flat on the sur- 
face 16 * . Electron micrographs of a wheat 
leaf and a corn leaf both reveal irregular 
acicular wax structures distributed on a 
flat wax surface much like tree stumps on 
a flat, clear cut area or a child's jacks 
scattered on a table 16 * *\ 

A. THEORY 

The mathematical model for polarisa- 
tion of light from a wheat canopy is de- 
veloped in two parts. First, the micro- 
scale situation, the polarization response 
and orientation requirements for a small 
specularly reflecting area Aa on a leaf 
(Fig. 3), is analyzed. Second, the micro- 
scale results are extended to the macro- 
scale level as measured by a field spec- 
trometer or satellite sensor (Fig. 4). 

The assumptions are: 

1. There exist on the wax surface of 
each shiny leaf small flat areas, Aa, 
which specularly reflect light. 

2. The wax layer is essentially clear 
and absorbs little light. This means 
that for the wax layer the complex index 
of refraction can be adequately approxi- 
mated by its real component, a reasonable 
supposition for the visible spectral re- 
gion where any light energy absorbed by 
the wax layer is then unavailable to the 
chloroplasts to promote photosynthesis. 
Limited evidence supports this assumptiorf 0 . 

3. Specular light reflection occurs 
principally at the air-wax boundary. Com- 
paratively negligible amounts of light are 
reflected specularly to an observer and 
from the boundaries between epidermal cell 
walls, cell membranes, and the various 
cutible layers. These boundaries have 
comparable indices of refraction and often 
appear rough in electron micrographs. 

4. The magnitude of polarized light from 
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sensor 

area A 


sources (moist soil, for example) other 
than sunlit leaves is insignificant. 

Micro Scale Response . On the micro 
scale level, Fig . 3 , sunlight is specular- 
ly reflected by one of the small leaf 
areas to an observer only if the area is 
properly oriented. The normal to the area 
must be in the plane and bisect the angle 
formed by two vectors, the illumination 
vector (e) directed from the leaf ar§a to 
the sun and the observation vector (V) di- 
rected from the leaf area to the observer. 
In other words the area Aa must be orient- 
ed such that the angle of incidence, y, 
equals the angle of reflectance and such 
that the vectors £,^n a (the unit vector 
normal to Aa) , and V are coplanar. Only 
then will specularly reflected light from 
Aa reach an obse-ver. These conditions 
form the keystor . of the polarization 
model and are satisfied by the equations, 

-+ -► -* -* 

E * n n -V 

— I a .id 

Y=cos =cos 

$\ |v| 

The solar irradiance incident on one 
small area, Aa, is a function of the angle 
of incidence. If the area Aa specularly 
reflects light to an observer, then the 
angle of incidence is uniauely determined 


=0 . 5COS* 1 


E-V 

w 


( 1 ) 



X toward 
sun 


Figure 3. Coordinate System. A 
small leaf area Aa specularly^ref lects 
sunlight toward an ojjsejver, V + if and 
only if the vectors E, n a and V are co- 
planar and the angles of incidence (y) 
and reflectance are equal. 



Figure 4. Canopy Response. A sensor 
measures the canopy response over a solid 
angle Aw. 

by the angles (0 S , 6 V , 4> v ) » as discussed 
above. For such an area the radiant flux 
incident on Aa is P s |e| cos y where y = 
y (S s , 6 V , i)> v ) and P s is the probability of 
finding in a small volume an area Aa il- 
luminated directly by the sun as opposed 
to being shaded by intervening foliage. 

The probability of finding in a small vol- 
ume an observable area Aa is symbolized by 
P v . In determining P s and P v , area Aa is 
assumed to be either illuminated or shaded, 
observable or not observable; adumbral 
effects are not considered 21 . 

The probabilities, P s and P v , are 
functions both of the (x,y,z) location of 
the leaf in the canopy and of directions 
of illumination and observations, respec- 
tively. The probabilities will be unity 
when only leaves are illuminated and ob- 
served. For example, P s and P v will ap- 
proach unity for the topmost leaves of a 
dense, preheaded wheat canopy if the ag- 
gregation of these leaves forms a layer 
one leaf thick at the extreme top of the 
canopy essentially impenetrable to direct 
illumination. The probabilities will be 
less than unity for more typical canopies 
with some soil and/or non-leaf foliage 
illuminated and observed and some leaves 
not illuminated and/or not observed. 

Even though the incident sunlight is 
not polarized, each small area, Aa, polar- 
izes the specular portion of the reflected 
sunlight provided the angle of incidence 
is neither 0 or 90 degrees. If Aa is 
smooth, the magnitude of the light that 
is specularly reflected and polarized by 
Aa is described mathematically by the 
Fresnel equations and Stokes vector and 
depends only upon the angle of incidence, 
Y, and the index of refraction of both 


100 


1980 /Machine Processing of Remotely Sensed Data Symposium 



the epicuticule wax layer and air". 

For a smooth surface the portions of 
the specular reflectance with electric vec- 
tor perpendicular, p*p I C Y ) # and parallel, 
Psp| | (y) to the plane of incidence are 
given by the Fresnel equations 22 . The 
first component of the Stokes vector, Sj, 
provides the magnitude of light reflected 
by a surface 22 . The term Sj is the sum of 
the specularly (Sg) and diffusely (S D ) re- 
flected light from the surface. The sec- 
ond component, Sq, is the portion of Sj 
linearly polarized by the surface. 


S I * 

S D * S S 

(2a) 

s s - 

<P -Pj 

(2b) 

S Q " 


(2c) 


[dimensionless] 


Most often Aa is not a perfectly 
smooth surface but instead supports small 
acicular structures. These structures 
diffusely scatter light which would other- 
wise be specularly reflected. To account 
for the reduced amount of light, specular- 
ly reflected by an area Aa which is not 
perfectly smooth, the Fresnel equations 
is modified by a factor, K. In qeneral K 
[dimensionless] is a function of the 
angles (6g,6 v ,$ v ), wavelength, and side 
(top or bottom) lateral position and di- 
rection on the leaf. However, evidence 
suqgests that for any one leaf the wax 
acicular structures in most cases are ho- 
mogeneously and isotropically dispersed 
across the leaf surface 18 . It is assumed 
here that K (0 S , 0 V , $ v , X ) is identical for 
all leaves and is not a function of lat- 
eral position and direction on the leaf 
surface. The value of K varies between 
zero and one. 

Define a probability density function 
f a (0,4>) for leaves such that the probabil- 
ity that any one of the leaf areas Aa is 
oriented within a solid angle Aw a about 
(^a'^a) is A(o a f a (e a ,<f a ) [dimensionless] . 

The units of f a =f a (x,y , z , 0 , if ) are [sr -1 ]. 
Because the area Aa must be correctly ori- 
ented to reflect light to an observer, the 
Jacobian provides 28 


area is the quadruple product of (1) the 
radiant flux incident on the area Aa, (2) 
the specular reflectance of Aa, (3) the 
probability that Aa is correctly oriented 
to specularly reflect light in a solid 
anqle Aw v about direction (6 v ,f v )i and (4) 
the probability that Aa is observed. 


( 1 ) ( 2 ) ( 3 ) ( 4 ) 

A4 G = P | £ | Aa cos> KS C Aw f (6 ,<f ) P (4) 

do S D a a a a V 

For a volume Vj, 4vjS = jjd^aS (5) 

all leaf area in Vj 


If the volume Vj is sufficiently small, 
then the probabilities P s and P v are essen- 
tially constants everywhere in Vj and will 
be denoted P s j and P v j. Letting A-j be the 
leaf area in volume Vj , f a j be the proba- 
bility density function of leaf area ori- 
entation in Vj, and usinq eq.(3) to sub- 
stitute for Aw a , eq.(5) becomes 

^Vj-Psj |E|AjKSgA v cos 2 e v f aj ( 0 a , 4> a ) P v j/ 4h 2 (6) 


Sur.minq the specular flux contributions of 
each volume Vj in the field of view of the 
sensor to find the specular portion of the 
flux measured, eq . 7, and the linearly po- 
larized portion of , eq.fi 

all Vj in 
field of view 


♦S-|E| KSjJ 

- f ai r si p VjAj 
= 1 

A v cos‘ 0 v /4h : 

source 

canopy 

sensor 

dependent 

dependent 

dependent 

term 

terms 

terms 

S’ 

^ ' 



all Vj in 
field of view 


♦q-|e| KS q £ fajP s i r v-i A i A v cos 2 e v /4h 2 
j*l J J 


(7) 


( 8 ) 


The percent linear polarization is propor- 
tional to the linearly polarized flux di- 
vided by the total flux, the sum of the 
diffuse and specular fluxes. 


Aw. = Aw /4cosy * A cos 2 0 /4h : cosv (3) 

« V V V 


1 00 %4q/ ( + 4g) 


( 9 ) 


Macro Scale Response . On the macro 
scale level. Fig. 4 , the radiant flux due 
to specularly reflected sunlight received 
from leaves in the field of view of a sen- 
sor is found by summinq the flux contrib- 
uted by each leaf area Aa in each volume 
Vj in the field of view. For a Aa the 
radiant flux specularly reflected into a 
solid anqle Aw v by a randomly selected 


To illustrate the properties of the polar- 
ization model, the response of several 
canopies will be examined. 

Example A: Sparse Wheat Canopy . I f 
the properties of the canopy are constants 
for those layers containing leaves, that 
is, if P s j=P s i=Ps' Pvj“Pvi*P v * A j=Ai = A, 
and faj = fai =f a' vv j» v i E field of view 
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Fiqure 5 . Preheaded Wheat Canopy 
Polarization Response. Prior to heading 
the response is zero at the anti-solar 
point, the "hot spot," and increases with 
increasing zenith view angle. 



Figure 6. Headed Wheat Canopy Polar- 
ization Response. After heading the re- 
sponse remains zero at the anti-solar 
point , is maximum .at intermediate zenith 
view angles, and approaches zero for near- 
horizontal view directions where heads and 
stems obstruct view of polarizinq flag 
leaves . 


and 

* s - | E | P 8 ( I .A I ) KS s f a P v A v .W2cos *0 v (10a) 


' — 



source 

canopy 

sensor 

. 


» , — . 


1 S'0*|p| r s< LAI ) KS 0 f a p v AwA v /2cos0 v (10b) 

where I.AI is leaf area index. The linear- 
ly polarized portion of the total scene 
radiance is*\ 

Lq-^q/AvA.)- |E| P s I.AIKS 0 f a P v /2cosO v (11) 

Example B: Pighe aded, Dense Wheat 

C anopy^ If the probability Psj p vj : l for 
V: in the topmost layer of the canopy and 
PgjP v j'0 for Vj in all lower layers, then 
tile •ts and Jq are proportional to the leaf 
area index only of the topmost layer. A 
winter wheat canopy measured just prior 
to heading might have the following char- 
acteristics: LAI ^2.0 for top layer con- 

taining flag leaves with a wax layer index 
of refraction * 1.5, f a ■ uniform 0.0796 
sr" 1 , K ■ 0.9, and P s = P v ■ 1.0 for top 
layer and P s = P v " 0.0 for all lower lay- 
ers. The linearly polarized portion of 
the canopy radiance (eg. 11) and the lin- 
early polarized flux measured by a sensor 
over such a canopy are shown in Fig. 5. 

The calculations are for a sensor with a 
field o' view. Aw, of 15° ; 0.216sr, en- 
trance optics of area 0.002m', and spec- 
tral band of 0.6-0. 7pm (red wavelengths). 
The solar insolation is assumed to be 
165.3W/m J in the 0.6-0. 7pm spectral band’’*. 

Exa mple C: Head ed, Dense Wheat Can - 

opy . The polarization response of a wheat 
canopy is expected to change significantly 
during the heading growth stage (Fig. 2). 
This is because the probabilities P s and 
F v of a headed wheat canopy, unlike those 
of a preheaded canopy, are pronounced 
functions of sun angle and view angle. 

The product r s P v may be estimated for a 
hypothetical headed canopy with LAI 2.0 
by applying linear regression techniques 
to data for a canopy with LAI *1.0 and 
scaling by a factor of 2. O'*. 

( 12 ) 

P s r v ».232exp(-1.17( ( 1 /cos0 v ) + ( l/cos0 s ) ) ) 


s . t . A^0 , then 

all Vj in 
field of view 

^ p s j p v-i A i 

j-1 1 

■ P s P v t total two sided leaf area in FOV] 
- P S P V [ 2 LAI ] [ h s Aw/cos 5 0u> v ] 


Equation 12, derived assuminq P s and P v to 
be independent, provides erroneous esti- 
mates of P S P V at angles near the canopy 
hot spot direction where the probability 
P(leaf observat ion | leaf illumination) 
approaches unity. Fig. 6 shows the linear 
polarization response for a source-canopy- 
sensor with the parameters of Example R 
except P S P V given by equation (12). 
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Examp Lr D: Percent Polarization 

Respo n se of Canopy^ If the magnitude of 
the specular flux is either small or larqe 
compared to the diffuse flux, then the 
percent polarization (eq. 9) reduces to 


* polarization 


1 0 0 % 't' q/ * ‘t’D s N (13a) 
100%'fg/*S'*D << ' ,, S (13b) 


Both Qq and ^g are proportional to the 
same agronomic factors, when t , D v< ' , *'S 
(which miqht occasionally be true in the 
blue and red spectral regions), then there 
is no aqronomic information in the term 
percent polarization because the aqronomic 
factors in the numerator and denominator 
of eq . 13b cancel. 


IV. DISCUSSION 

The model shows the magnitude of the 
polarization response of a plant stand is 
related to the solar insolation and the 
characteristics of the canopy and the sen- 
sor. The response depends on the optical 
and qeometric properties of the portion of 
the canopy in the instrument field of view. 

The calculations show that a sensor 
would measure zero linearly polarized 
liqht at the anti-solar point, the canopy 
"hot-spot," where Sq and the angle of in- 
cidence of the sunlight are both zero. 

The sensor would measure the maximum 
amount of polarized liqht in the solar 
azimuth direction (Fiq. 5) , provided the 
small areas Aa are randomly oriented in 
azimuth and zenith directions. Otherwise 
the direction of maximum polarization may 
be shifted, as might occur when a strong 
wind preferentially orients the flag 
leaves of wheat downwind. 

The theory shows that when the ap- 
proximation (eq. 13b) is valid--when the 
specular flux is much, much qreater than 
the diffuse flux — the percent linear po- 
larization is not directly related to the 
canopy agronomic properties. This might 
occur in the chlorophyll absorption re- 
gion in the red portion of the spectrum 
viewinq at larqe zenith anqles toward the 
sun azimuth angle. However, even though 
the percent linear polarization in certain 
circumstances may contain limited infor- 
mation related to canopy agronomic fac- 
tors, the magnitude (eq. 8) of the linear- 
ly polarized flux is always directly re- 
lated to the canopy aqronomic properties. 
Thus, the model provides a theoretical 
basis for the same, but empirically based 
result noticed by Egan’. 

The aqronomic variables LAI, f a , K, 
Sq, P s , and P v in the equations (eqs. 8 


and 10b) are functions of one or more en- 
vironmental and/or physiological variables. 
The leaf area index (LAI) is a function of 
many variables including species and cul- 
tivar, weather, and growth stage. The 
probability density function f a (0 a »$ a ) f° r 
the orientation of the leaf areas Aa is a 
function of wind strenqth and direction, 
catastrophes (such as hail damage) , crop 
vigor (moisture stress causes corn leaves 
to roll, cotton leaves to droop), and 
qrowth staqe (the shapes of healthy green 
leaves and senescent leaves are not the 
same) . The optical properties of the 
leaves, K, Sq, and Sj# are functions of 
species and almost certainly cultivar, di- 
sease (plant pathogens often alter or de- 
stroy the wax layer) , pubescence (the 
hairs scatter light which would otherwise 
be specularly reflected) , material on the 
leaves (dust, pollen, water droplets), and 
wavelength (the Fresnel equations are func- 
tions of the cuticle wax index of refrac- 
tion, which changes in regular fashion 
with wavelength) . The terms Sq and Sg are 
functions of the angle of incidence (view 
and illumination direction) of the sun- 
light. The geometrical properties of the 
canopy, P s and P v , are functions of the 
anqles of illumination and observation 
(lower leaves have lower probabilities of 
beinq illuminated and observed at large 
sun and view zenith angles) and growth 
stage (wheat heads partially block illumi- 
nation and observation of flag leaves; the 
projected area of leaves changes with se- 
nescence) . Several of these functional 
relationships will be discussed further 
in the following paragraphs. 

Liqht polarized by a moist or wet 
soil surface is a part of the canopy po- 
larization response not considered in the 
theory. Visual observations suqgest that 
except for wet soil surfaces, the amount 
of liqht polarized by the soil is insig- 
nificant compared to the amount of light 
polarized by foliage. Neglecting emergent 
and sparse canopies, the soil generally 
has a very low probability of being both 
illuminated and observed. Thus, the the- 
ory presented predicts the polarization 
response of canopies on dry soils and/or 
with sufficient foliaqe to obscure the 
soil . 


Sunlight tends to be specularly re- 
flected and polarized by leaves in the 
upper portion of plant stands. The proba- 
bility that a leaf is both visible and 
directly illuminated by sunliqht is often 
a pronounced function of the (x,y,z) lo- 
cation of the leaf in the stand and of 
the view and illumination directions. 

The probability tends to decrease rapidly 
with increasing depth into the canopy. 

This means that leaves in the lower portions 
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of a plant stand will little affect the 
canopy polarization response. Thus, lower 
leaf senescence or a disease condition lo- 
calized to the lower leaves may not be de- 
tectable using polarization measurements. 

The probability density function, 
fa^an^a), for the orientation of the leaf 
areas, Aa, can be calculated (from eq. 10b) 
for the population of observable, specu- 
larly reflecting leaves using polarization 
measurements sampled from the hemisphere 
of all possible canopy view directions. 
Such a density function is needed as input 
data to certain canopy radiation models 
which are used to examine the utility of 
and information in canopy reflectance 
measurements*’ 7 . To obtain the polarization 
data needed to calculate f a , Horvath pro- 
posed a field apparatus consisting of a 
linearly polarized light source co-located 
with a sensor with a polarization analy- 
zer*’*. Due to the inherent randomness of 
the leaf structure, the diffuse portion 
of the reflected light will tend not to be 
polarized". The specular portion of the 
reflected light will be polarized, not 
because of reflection at the leaf cuticle 
wax layer (0° angle of incidence) but in- 
stead because the liaht source is polar- 
ized. To compute f a the canopy polariza- 
tion response is measured with the analy- 
zer oriented in two directions, parallel 
and perpendicular to the polarization di- 
rection of the light source beam. It is 
easily shown that the specularly reflected 
flux (eq. 10a) is the difference of these 
two measurements; hence, in a particular 
direction 

f a = |E|LAI KS s P s P v AwA v /21' s cosO v 

where Sg = ( (n-1) / (n+1 ) ) *’ for normal inci- 
dence, and P S =P V . A practical limitation 
to the approach exists. In general values 
for leaf area index (LAI), the factor K, 
and the index of refraction (n) , are not 
known; however, properly normal izinq to 
unity the integral with (0 v ,<t v ) of the ini- 
tial estimate of f a obviates the need to 
know these terms. But more importantly 
the need to know or estimate Pg and P v , 
usually unknown functions of view direc- 
tions, cannot be circumvented if f a is to 
be calculated. The term P s p yfa includes 
all the canopy dependent variation due to 
view angle (assuming the factor K is a 
constant) and is always calculable. 

Depending on its direction and 
strength, the wind is capable of reorient- 
ing the leaves of a canopy and thereby 
changing the probability density function 
of leaf area, f a . The resultant varia- 
tions from day to day in the polarization 
response of a field will tend to compli- 
cate interpretation of polarization data 


because these variations represent noise 
(unless, of course, the probability densi- 
ty function of the orientation of the 
leaves is to be calculated for each day) . 
The size of these day to day variations, 
if sufficiently large, might preclude a 
naive analysis which neglects wind effects; 
conceivably, wind induced variations in f a 
might render an agronomic interpretation 
of polarization data impossible. There 
remains the hope, however, that for data 
taken at one time the wind will affect 
similarly all the fields in a region con- 
taining a particular crop species and cul- 
tivar. For such data comparisons between 
fields of a specific crop variety might 
remain valid. Therefore, an important 
question which should be addressed empiri- 
cally is the following; How uniformly 
does wind affect polarization data acquired 
over a region? 

Visual observations suggest that blue 
skylight incident on the canopy affects 
minimally the magnitude of the canopy po- 
larization response. By the same process 
discussed in the theory for sunlight, the 
shiny leaves of a canopy polarize the 
specularly reflected skylight, a spectrally 
varying light source already polarized 
according to observer view direction. The 
magnitude of the skylight and its effect 
on the canopy polarization response is 
qreatest in the blue spectral region and 
decreases into the near infrared. Atmos- 
pheric haze, which decreases the solar 
insolation on leaves, noticeably decreases 
the specular and polarized light from 
leaves . 

The efficacy and feasibility of a 
satellite sensor measuring the linear po- 
larization of a scene through the earth's 
atmosphere has not been considered. Sen- 
sor design must consider the path radiance 
of the atmosphere, a source of linearly 
polarized light potentially capable of 
altering or masking the amount of linearly 
polarized light received from the scene. 
Even if the polarized portion of the path 
radiance is excluded, analysis of field 
spect-al radiometer data suggests that for 
a satellite polarization sensor the signal 
flux must be increased (the spectral/spa- 
tial resolution of a Landsat-type satel- 
lite must be degraded probably by a factor 
greater than 10:1) and/or noise power de- 
creased to obtain a signal to noise ratio 
approximately equivalent to that of Land- 
sat in the red spectral channel. Including 
atmospheric effects in the analysis would 
potentially indicate a practical value of 
the spatial/spectral resolution of a sat- 
ellite sensor. 

The information in canopy polarization 
data, when obtained from satellite sensors, 
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Figure 7. Polarization Response with Heading. The three photographs, taken at one 
week intervals before and after heading, demonstrate that the amount of specularly re- 
flected- and therefore polarized-sunlight decreases significantly with heading. 


probably will be used in conjunction with 
other remotely sensed data and will be 
extracted by analysis of frequent, synop- 
tic data sets, by using the temporal and 
spatial information to make relative com- 
parisons between the fields in the data 
for one date and between the dates for 
one field. One polarization measurement 
of one field for one date probably will 
have little value unless it is compared 
to polarization data for that field and 
other fields for that date and other dates. 
This is because it is unrealistic from the 
model to expect that canopy polarization 
data will be calibrated in an absolute 
sense to discriminate a particular crop 
or to correlate uniquely to a particular 
agronomic variable. Frequent, synoptic 
polarization data from a satellite sensor 
potentially aid in assessing crop vigor 
and growth stage and in determining areas 
of hail damage and pestilence, all poten- 
tially possible from comparisons between 
field and across dates. Perhaps daily 
satellite coverage is feasible using a low 
spatial and spectral resolution sensor in 
a geosynchronous orbit. 

The canopy polarization response de- 
scribed by the model is a function of 
wavelength only because the index of re- 
fraction of the cuticle wax layer is a 
function of wavelength. From the physics 
of the optical properties of materials it 
is expected that the index of refraction 
of the wax layer will gradually and mono- 
tonically increase with decreasing visible 
wavelength, displaying no perturbations 
or "fine structure" with wavelength 14 . 
However, the model indicates the percent 
linear polarization of a healthy green 
canopy will be large in the blue and red 
spectral regions, small in the green, and 
even smaller in the near infrared region 
away from any absorption bands. This is 
because the total canopy flux, the normal- 


ization factor used when computing percent 
polarization, exhibits a green vegetation 
response 2 \ 

From the model there appears little 
need to measure the canopy spectral polar- 
ization response with high wavelength re- 
solution in the visible spectral region; 
a polarization sensor covering the entire 
visible region or a large portion of it 
might suffice. Conversely, in the infra- 
red spectral region the cuticle wax layer 
may absorb in narrow spectral regions de- 
fined by the structural properties of the 
constituent waxes of the layer, by the 
resonant frequencies of the translational 
and rotational vibration modes of mole- 
cules of the layer 14 . If absorption bands 
exist, high resolution spectral polariza- 
tion data may possibly provide information 
concerning the properties of the cuticle 
wax layer, properties relatable to crop 
species and light regime 18 . 

The linear polarization model may be 
extended to include the elliptical polar- 
ization response of the canopy. Evidence 
exists that the cuticle wax of some species 
is birefringent and therefore potentially 
capable of elliptically polarizing specu- 
larly reflected light. However, this 
evidence does not suggest that the struc- 
ture of the wax layer is sufficiently 
organized to elliptically polarize light 
specularly reflected from a significant 
amount of leaf area. Egan argues that 
the amount of elliptically polarized light 
from a vegetation scene should be negli- 
gible because of the inherent randomness 
of the properties of the vegetation 9 

The connection between leaf polariza- 
tion measurements and leaf moisture con- 
tent, noted by Egan, is supported by a 
morphological model for the structural 
changes which occur in a leaf undergoing 
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dehydration”. When the leaf water content 
and leaf thickness decrease and the leaf 
cells dehydrate and collapse, the surface 
roughness of the wax layer increases. The 
leaf surface appears rough in thin sections 
of senescent leaves and of leaves under 
moisture stress' 0 . As the surface rough- 
ness increases, the specular portion of 
the light reflected by the leaf decreases 
because there are fewer areas Aa which are 
similarly oriented. The amount of linear- 
ly polarized light reflected by the canopy 
decreases in company with the decrease of 
specularly reflected light. These argu- 
ments suggest that the canopy polarization 
response should decrease with decreasing 
leaf water content in the canopy and 
therefore serve as an indicator of canopy 
moisture stress. Visual evidence supports 
this hypothesis. Leaves under moisture 
stress often appear less shiny than fully 
turgid leaves. Dry, senescent leaves 
often have a matte surface finish. 

Detection of the date of heading of 
a wheat canopy (Fig. 2), information which 
is needed for use with phenologically 
based models to predict the ultimate grain 
yield of the crop, might be feasible using 
satellite polarization measurements (Fig. 

7) . The eventual weight of grain produced 
by each wheat plant is largely determined 
by the condition of the flag ieaf, its 
size and vigor, and by the weather regime 
endured by the plant following heading 
when the grain head begins to fill 30 . Know- 
ledge of the date of heading permits a 
better estimate of the post-heading weather 
for the crop. Prior to heading the top- 
most foliage on the wheat plant is the 
flag leaf, easily the most visible and 
illuminated canopy component (Fig. 2). 
Following heading, wheat heads are the 
topmost foliage and partially obscure the 
flag leaves to both sunlight and observa- 
tion, changing the values of both P s and 
P v . Figures 5 and 6 show that the magni- 
tude of the polarized light, which depends 
directly upon the specular reflections 
from flag leaves, will decrease by a fac- 
tor of 60 for 0 S =3O and 0y=O for the two 
hypothetical canopies during heading as 
the leaves are increasingly obscured to 
both illumination and observation. The 
obscuration of the flag leaves is enhanced 
at off nadir observation angles directed 
toward the solar azimuth (Fig. 6) . Poten- 
tially both the condition of the flag 
leaves and the date of heading of a crop 
might be monitored using polarization 
measurements obtained from a satellite 
sensor with both on and off nadir viewing 
capability. Such a view capability has 
been proposed for the MRS sensor. 

Applicability of the polarization 
model should extend to many species because 


shiny leaves which specularly reflect sun- 
light are ubiquitous, unconfined by geog- 
raphy or climate. Other plants besides 
wheat, sorqhum, and corn with specularly 
reflecting leaves include coffee, sudan 
grass, banana, orange, sugarcane, and many 
forest species. Schieferstein and Loomis 
found epicuticular wax deposits on about 
half of the plant species they tested 1 ’. 
However, the mere presence of a cuticle 
wax layer does not guarantee that a leaf 
will specularly reflect and polarize a 
significant portion of the incident light; 
the leaf must also appear shiny. Fibrillar 
light scattering significantly diminishes 
the polarization response of pubescent 
soybean leaves. And the surface of the 
wax layer of some species is insufficient- 
ly smooth to specularly reflect light 18 . 


V. CONCLUSIONS 

This paper discusses a model for the 
amount of linearly polarized light reflect- 
ed by* the shiny leaves of such crops as 
wheat, corn, and sorghum, each a grain of 
major economic importance to the world. 

The model is based upon the morphological 
and phenological characteristics of the 
canopy and upon the Fresnel equations 
which describe the light reflection pro- 
cess at the smooth boundary separating two 
dielectrics . 

The theory demonstrates that, poten- 
tially, measurements of the linearly po- 
larized light from a crop canopy may be 
used as an additional feature to discrim- 
inate between crops such as wheat and 
barley, two crops so spectrally similar 
that they are misclassif ied with unaccep- 
table frequency. Examination of the model 
suggests that, potentially, satellite po- 
larization measurements may be used to 
monitor crop development stage, leaf water 
content, leaf area index, hail damage, and 
certain plant diseases. Such information 
is needed for use with models which pre- 
dict crop grain yield. 

The model adds to our understanding 
of the potential information content of 
scene polarization measurements. The in- 
formation content of these measurements 
has not been extensively investigated and 
needs to be understood to evaluate the 
potential usefulness of the proposed po- 
larization sensor for the satellite borne 
Multispectral Resource Sampler. The effi- 
cacy of a satellite sensor measuring the 
linear polarization of a scene through the 
atmosphere remains to be determined. 
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